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Abstract 

The Stable Marriage Problem (SMP) is a well-known matching problem first introduced and 
solved by Gale and Shapley [7]. Several variants and extensions to this problem have since been 
investigated to cover a wider set of applications. Each time a new variant is considered, however, 
a new algorithm needs to be developed and implemented. As an alternative, in this paper we 
propose an encoding of the SMP using Answer Set Programming (ASP). Our encoding can easily 
be extended and adapted to the needs of specific applications. As an illustration we show how 
stable matchings can be found when individuals may designate unacceptable partners and ties 
between preferences are allowed. Subsequently, we show how our ASP based encoding naturally 
allows us to select specific stable matchings which are optimal according to a given criterion. Each 
time, we can rely on generic and efficient off-the-shelf answer set solvers to find (optimal) stable 
matchings. 

Keywords: Answer Set Programming, Logic Rules, Stable Marriage Problem, Optimal Stable Matchings. 

1 Introduction 

The Stable Marriage Problem (SMP) is a matching problem first introduced and solved by Gale and 
Shapley [?]■ Starting from (i) a set of n men and n women, (ii) for each man a ranking of the 
women as preferred partners, and (iii) for each woman a ranking of the men as preferred partners, 
the SMP searches for a set of n couples (marriages) such that there are no man and woman who are 
in different marriages but both prefer each other to their actual partners. Such a man and woman 
are called a blocking pair and a matching without blocking pairs forms a stable set of marriages. Due 
to its practical relevance, countless variants on the SMP have been investigated, making the problem 
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assumptions more applicable to a wider range of applications, such as kidney-exchange [12] and the 
hospital-resident problem [TB]. Recently Roth and Shapley won the Nobel Prize for Economics for 
their theory of stable allocations and the practice of market design, work that has directly resulted 
from an application of the SMP. 

In the literature, typically each time a new variant on the SMP is considered, a new algorithm is 
developed (see e.g. [Tni[T31[^). In this paper, we propose to use Answer Set Programming (ASP) 
as a general vehicle for modeling a large class of extensions and variations of the SMP. We show how 
an ASP encoding allows us to express in a natural way ties in the preferences of men and women, as 
well as unacceptability constraints (where certain people prefer to remain single over being coupled 
to undesirable partners). Furthermore, we illustrate how we can use our ASP encoding to find stable 
matchings that are optimal according to a certain criterion. Although the SMP has been widely 
investigated, and efficient approximation or exact algorithms are available for several of its variants 
(see e.g. [IS]), to the best of our knowledge, our encoding offers the first exact implementation to find 
sex-equal, minimum regret, egalitarian or maximum cardinality stable sets for SMP instances with 
unacceptability and ties. 

The paper is structured as follows. In Section[2]we give some background about the SMP and ASP. 
We introduce our encoding of the SMP with ASP and prove its correctness in the third section. In 
Section 131 we extend our encoding enabling it to find optimal stable sets. We explore several notions 
of optimality for stable matchings and show how optimal stable matchings can be found by solving 
the corresponding disjunctive ASP program. Finally we draw our conclusions. 

2 Background 

2.1 The Stable Marriage Problem 

To solve the standard SMP, Gale and Shapley ^ constructed an iterative algorithm — known as the 
Gale-Shapley algorithm, G-S algorithm or deferred-acceptance algorithm — to compute a particular 
solution of an SMP instance. The algorithm works as follows: in round 1 every man proposes to his 
first choice of all women. A woman, when being proposed, then rejects all men but her first choice 
among the subset of men who proposed to her. That first choice becomes her temporary husband. 
In the next rounds, all rejected men propose to their first choice of the subset of women by whom 
they were not rejected yet, regardless of whether this woman already has a temporary husband. Each 
woman, when being proposed, then rejects all men but her first choice among the subset of men 
who just proposed to her and her temporary mate. This process continues until all women have a 
husband. This point, when everyone has a partner, is always reached after a polynomial number of 
steps and the corresponding set of marriages is stable [7] . It should be noted, however, that only one 
of the potentially exponentially many stable matchings is found in this way. We formally define the 
SMP and introduce two variants that will be considered in this paper. We denote a set of men as 
M = {toi, . . . , m„} and a set of women W = {wi, . . . , Wp}, with n = p for the classical SMP. A set of 
marriages is a set of man-woman pairs such that each man and each woman occurs in just one pair. 

Definition 1 (Classical SMP) 

An instance of the classical SMP is a pair {Sm, Sw), with Sm — ■ ■ ■ , o'm} and Sw — Wwj • ■ • : ^w} 
sets of permutations of the integers 1, . . . , n. The permutations and cr^ are the preferences of man 
rui and woman Wi respectively. If a\j{j) — k, we say that woman Wk is the j*^ most preferred woman 



for man mi, and similarly for cr^(j) — k. Man m and woman w form a blocking pair in a set of 
marriages S if m prefers w to his partner in S and w prefers m to her partner in S. A solution of 
an instance is a stable set of marriages, i.e. a set of marriages without blocking pairs. 

A first variant of the classical SMP allows men and women to point out unacceptable partners 
by not including them in their preference list. The number of men n can differ from the number of 
women p since men and women can remain single. A set of marriages is a set of singles (i.e. persons 
paired to themselves) and man-woman pairs such that every man and woman occurs in just one pair. 

Definition 2 (SMP with unacceptability) 

An instance of the SMP with unacceptability is a pair {Sm,Sw), Sm — Wmt--t'^m}, o,i^d Sw — 
{crjy, . . ., o'yy}, with cuch a permutation of a subset of {1, . . . ,p} and each cr^ a permutation of 
a subset of {1, . . . ,n}. If cr\.i{j) — k, woman Wk is the j*^ most preferred woman for man mi, and 
similarly for crfy(j) = k. If there is no I such that (t\,j{1) — j, woman Wj is an unacceptable partner 
for man mi, and similarly for no I such that (j\y{1) ^ j. A person x forms a blocking individual in a 
set of marriages S if x prefers being single to being paired with his or her partner in S. A solution of 
an instance is a stable set of marriages, i.e. a set of marriages without blocking pairs or individuals. 

The length of the permutation a\,j is denoted as \(J%i\. A stable matching for an SMP instance 
with unacceptability always exists and can be found in polynomial time [22] by a slightly modified 
G-S algorithm. 

Example 1 

Suppose M = {mi, 7712, mg}, W = {^1,^2,^3,^4}, Sm = {crli = (4, 1, 3), ct^.j (3,2), af j = (1,3)} 
and Sw = Ww ~ (1,3),(t^ = (2),(7^ = (3,2),(T^ — (2,1)}. Hence woman wi prefers man mi to 
man m^ while man m2 is unacceptable. In this setting, there is exactly one stable set of marriages |22j : 
{(mi, u;4), (m2, W3), (m3, wi), (i(;2, W2)}. Thus woman ^2 stays single. 

The second variant of the SMP allows unacceptability and ties, i.e. the preferences do not have to 
be strict. For this variant there are several ways to define stability, but we will use the notion of weak 
stability |11) . 

Definition 3 (SMP with unacceptability and ties) 

An instance of the SMP with unacceptability and ties is a pair (^S ]\j , S\y^ , S]\j — {fTj|j,...,fj^^} and 
Sw = Ww^ • ■ • ' '^w}- ^^''^ every i £ {1, . . . , n}, a\.i is a list of disjoint subsets of {1, ... ,p}. Sym- 
metrically is a list of disjoint subsets 0/ {1, . . . , n} for every i £ {1, . . . ,p}. We call a\j and a\y 
the preferences of man mi and woman Wi respectively. If k € (^^^{j), woman Wk is in man mi 's j*'' 
most preferred group of women. All the women in that group are equally preferred by mi. The case 
k £ cr|^(j) is similar. If there is no I such that j £ a\j{l), woman wj is an unacceptable partner for 
man mi, and similar for no I such that j £ a\y{l). For every k in the se£] cr\^(|CT\^|), man mi equally 
prefers staying single to being paired to woman Wk, and symmetrically for the preferences of a woman 
Wi . This is the only set in a]^ that might be empty, and similar for . Man m and woman w form 
a blocking pair in a set of marriages S if m strictly prefers w to his partner in S and w strictly prefers 
m to her partner in S . A blocking individual in S is a person who stricly prefers being single to being 
paired to his partner in S. A solution of an instance is a weakly stable set of marriages, i.e. a set of 
marriages without blocking pairs or individuals. 

^|<t]^| denotes the length of the hst cr^y^. 



A weakly stable matching always exists for an instance of the SMP with unacceptability and ties 
and it can be found in polynomial time by arbitrarily breaking the ties |14| . However, as opposed to 
the previous variant, the number of matched persons is no longer constant for every stable set in this 
variant. Note that the setting of Definition |3] generalizes the setting of Definition [2l which generalizes 
the setting of Definition [TJ We introduce the notations 



acceptable 




— preferred\.j — neutral\j 



Furthermore unacceptable\j = {l,...,p} \ acceptable\.j. We define the ordening on {wj \j G 
acceptable\j} U {rrii} as x y iff mi prefers person x at least as much as person y. The strict 
ordening <™j is defined in the obvious way and analogous notations are used for a^y. 

Example 2 

Suppose M = {mi, 7712}, W = {wi, W2, W3, W4} and Sm = Wlj = ({1,3}, {4}), ajj = ({2,3},{})}. 
Hence man mi prefers women wi and lus to woman w^. There is a tie between woman wi and ws 
as well as between woman and staying single. Woman W2 is unacceptable for man mi. Man m2 
prefers woman W2 and W3 to staying single, but finds wi and unacceptable. It holds that wi <^^/ 
mi, i.e. mi prefers marrying wi over staying single, acceptable\.j — {1,3,4}, preferred\j — {1,3}, 
neutral\j — {4} and unacceptable]^ = {2}. 



2.2 Answer Set Programming 

Answer Set Programming or ASP is a form of declarative programming [2] . Its transparence, elegance 
and ability to deal with -complete problems make it an attractive method for solving combinatorial 
search and optimization problems. An ASP program is a finite collection of first-order rules 

Ai V . . . V Afc ^ Bi, . . . , B™, notCi,..., not Cn 

with Ai, . . . , Ak, i?i, . . . , B,n, Ci, . . . , Cn predicates. The semantics are defined by the ground version 
of the program, consisting of all ground instantiations of the rules w.r.t. the constants that appear in 
it (see e.g. [2j for a good overview). This grounded program is a propositional ASP program. The 
building blocks of these programs are atoms, literals and rules. The most elementary are atoms, which 
are propositional variables that can be true or false. A literal is an atom or a negated atom. Beside 
strong negation, ASP uses a special kind of negation, namely negation-as-failure (naf), denoted with 
'not'. For a literal a we call ^nota^ the naf-literal associated with a. The extended literals consist of 
all literals and their associated naf-literals. A disjunctive rule has the following form 

ai V . . . V afc ^ 61, . . . , bm,notci, ...,notCn 

where ai, at, 61, ...,&,„, ci, c„ are literals from a fixed set C, determined by a fixed set A of 
atoms. We call oi V . . . V Ofc the head of the rule while the set of extended literals 61 , . . . , fe„i , not ci , . . ., 
notCn is called the body. The rule above intuitively encodes that ai, 02, ... or is true when we 
have evidence that bi, . . . ,bm are true and we have no evidence that at least one of ci, . . . , c„ are true. 
When a rule has an empty body, we call it a fact; when the head is empty, we speak of a constraint. A 
rule without occurrences of not is called a simple disjunctive rule. A simple disjunctive ASP program 
is a finite collection of simple disjunctive rules and similarly a disjunctive ASP program is a finite 
collection of disjunctive rules. If each rule head consists of at most one literal, we speak of a normal 
ASP program. 



We define an interpretation / of a disjunctive ASP program "P as a subset of C An interpretation / 
satisfies a simple disjunctive rule oiV. . .Vofc <— &i, . . . , 6™ when ai G /V. . .Vat- £ / or {bi, . . . , &„} ^ /. 
An interpretation which satisfies all rules of a simple disjunctive program is called a model of that 
program. An interpretation / is an answer set of a simple disjunctive program V iff it is a minimal 
model of V, i.e. no strict subset of / is a model of 'P [9]. The reduct of a disjunctive ASP program 
V w.r.t. an interpretation / is defined as the simple disjunctive ASP program = {ai V ... V 
61, . . . , 6m I (fli V . . . V Ofc ^ fei, . . . , 6m, not Ci, . . . , not c„) 6 P, {ci, . . . , c„} fl / = 0}. An interpretation 
/ of a disjunctive ASP program V is an answer set of P iff / is an answer set of . 

Example 3 

Let V be the ASP program with the following 4 rules: 

man{john) <— , person{john) <— , person{fiona) <— 
woman{X) V child{X) ^ person{X),notman{X) 

The last rule is grounded to 2 rules in which X is resp. replaced by john and by fiona. We check 
that the interpretation / — {man(john),'woman{ fiona), person{john), per son{ fiona)} is an answer 
set of the ground version of V by computing the reduct. The grounded rule with X ~ john is deleted 
since manijohn) is in /. The reduct is; 

man{john) 4—, person{john) , per son{f iona) <— 

woman(fiona) V child{fiona) person{fiona) 

The first 3 rules are facts, hence their heads will be in any answer set. The fourth rule encodes that 
any person who is not a man, is a woman or child. It is clear that / is a minimal model of this simple 
program, so / is an answer set of V. By replacing woman{ fiona) by child{f iona) in /, another answer 
set is obtained. 

To automatically compute the answer sets of the programs in this paper, we have used the ASP 
solver DL\i|, due to its ability to handle predicates, disjunction and numeric values (with some built-in 
aggregate functions). The numeric values are only used for grounding. 

3 Modeling the Stable Marriage Problem in ASP 

In this section we model variations and generalizations of the SMP with ASP. A few proposals of 
using nonmonotonic reasoning for modeling the SMP have already been described in the literature. 
For instance in [19] a specific variant of the SMP is mentioned (in which boys each know a subset of 
a set of girls and want to be matched to a girl they know) and in [4 an abductive program is used to 
find a stable set of marriages in which two fixed persons are paired, with strict, complete preference 
lists. To the best of our knowledge, beyond a few specific examples, no comprehensive study has been 
made of using ASP or related paradigms in this context. In particular, the generality of our ASP 
framework for weakly stable sets of SMP instances with unacceptablity and/or ties is a significant 
advantage. The expression accept{m, w) denotes that a man m and a woman w accept each other 
as partners. The predicate manpropose{m, w) expresses that man m is willing to propose to woman 
w and analogously 'womanpropose{m,w) expresses that woman w is willing to propose to man m. 
Inspired by the Gale-Shapley algorithm, we look for an ASP formalisation to find the stable sets. 



^Available from www.dlvsystems.com 



Definition 4 (ASP program induced by SMP with unacc. and ties) 

The ASP program V induced by an instance {{clj, . . . , crjf/}, {crly, ■ ■ ■ , o'h/}) of the classical SMP with 
unacceptability and ties is the program containing for every i G {1, . . . , n}, j £ {1, . . . ,p} the following 
rules: 

accept(mi, Wj) manpropose{mi, Wj), womanpropose{mi, Wj) (1) 
accept{mi,mi) ^ {not accept{mi, w^) \ k G acceptable\.f\ (2) 
accept{wj,Wj) {not accept{mk, Wj) \ k G acceptable-!^} (3) 
and for every i G {!,..., n\ , j G acceptable'j^j : 

manpropose{mi, Wj) <— {not accept{mi, x) \ x Wj and Wj ^ x} (4) 
and for every j G {1, . . . ,p}, i G acceptablcy^ : 

womanpropose{mi,Wj) ^ {not accept{x,Wj) \ x <^ and ^ x} (5) 

Intuitively ([1} means that a man and woman accept each other as partners if they propose to each 
other. Due to a man accepts himself as a partner (i.e. stays single) if no woman in his preference 
list is prepared to propose to him. Rule Q states that a man proposes to a woman if he is not paired 
to a more or equally preferred woman. For j G neutral\j the body of Q contains notaccept(mi,mi). 
No explicite rules are stated about the number of persons someone can propose to or accept but 
Proposition [T] implies that this is unnecessary. 

We illustrate the induced ASP program with an example. 

Example 4 

Consider the following instance {Sm,Sw) of the SMP with unacceptability and ties. Let M = 
{mi, 7712} and W — {wi, W2,W3}. Furthermore: 





= ({1},{2,3},{}) 




= ({2},{1}) 




= ({1,2},{}) 


2 


= ({!},{}) 


^3 


= ({2}, {!},{}) 



The ASP program induced by this SMP instance is: 
man{mi) man(m2) ^ 

woman{wi) ,woman{w2) ■'r-, woman{w^) -ir- 

accept{X, Y) -S— manpropose{X , Y), womanpropose{X , Y), man{X), woman(Y) 
manpropose{mi,wi) -S— 

manpropose{mi,W2) notaccept{mi,wi),notaccept{mi,W3) 
manpropose{mi, W3) <~ notaccept{mi,wi),notaccept(rni,W2) 

accept{mi, mi) ^ not accept{mi,wi),not accept{m,i, W2), not accept {mi,w 3) 
manpropose{m2,W2) 

manpropose{m2,wi) ^ notaccept{m2,W2),notaccept(rn2,rn2) 
accept {m2,m2) ^ notaccept{m2,W2),notaccept(rn2,Wi) 



womanproposeijnijWi) -f- notaccept{m2,wi) 
womanpropose{m2,'Wi) ^ not accept{mi, wi) 

accept {w I, wi) ^ not accept{mi,wi),not accept{m2,wi) 
womanpropose{mi,W2) 

accept('W2,W2) not accept(mi, W2) 
womanpropose(m2,w^) -h- 
womanpropose{mi,W3) <— not accept{m2, w^) 

accept{w^,w^) not accept{mi, w^) , not accept{m2, w^) 

Notice that we use the facts man and woman to capture ah the rules of the form ([T]) at once. If we 
run this program in DLV, we get three answer sets containing respectively: 

• {accept{mi, w^) , accept{m2, wi) , accept{w2, 'W2)} , 

• {accept{mi, W2) , accept{m2, Wi) , accept{w3, w^)}, 

• {accept{mi, Wi),accept{m2,m2), accept{w2, W2), accept{w3, W3)}. 

These correspond to the three weakly stable set of marriages of this SMP instance, namely {(mi, W3), 
(m2,wi), (^2,^2)}, {(mi, W2), (m2,wi), (wa^ws)} and {(mi,wi), (m2,m2), (^2,^2), (103,^3)}. 

Proposition 1 

Let (SmtSw) be an instance of the SMP with unacceptability and ties and let V be the corresponding 
ASP program. If I is an answer set of V , then a weakly stable matching for {Sm,Sw) is given by 
{{x,y) I accept{x,y) £ /}. 

Proof. Let {Sm, Sw) and V be as described in the proposition. Because of the symmetry between the 
men and the women we restrict ourselves to the male case when possible. We prove this proposition 
in 4 steps. 

1. For every i G {1, . . . ,n}, every j G {1, . . . ,p} and for every answer set I of V , it holds that 
accept(mi,Wj) 6 / implies that j £ acceptable\j and i G acceptablc'l^ . 

This can be proved by contradiction. We will prove that for every man m^ and every j G 
unacceptable^]^, accept{mi,'Wj) is in no answer set / of the induced ASP program V. For 
accept{mi,Wj) to be in an answer set /, the reduct must contain some rule with this literal in 
the head and a true body. The only rule which can make this happen is the one of the form ([!]), 
implying that manpropose{mi,Wj) should be in /. But since j is not in acceptable\.j there is no 
rule with manpropose{mi,Wj) in the head and so manpropose{mi, Wj) can never be in /. 

2. For every answer set I of V and every man mi, there exists at most one woman Wj such that 
accept{mi,Wj) G /. Similarly, for every woman Wj there exists at most one man mt such that 
accept{mi,Wj) G /. Moreover, if accept{mi,mi) G / then accept(mi,Wj) ^ / for any Wj, and 
likewise when accept(wj, Wj) G / then accept{mi, Wj) ^ I for any mi. 

This can be proved by contradiction. Suppose first that there is an answer set I oi V that 
contains accept (mi,Wj) and accept{mi,Wji) for some man m^ and two different women Wj and 
Wj'. The first step implies that j and j' are elements of acceptable\j. Either man m^ prefers 
woman wj to woman Wji (wj <^j' Wji) or the other way around {wji Wj) or man m^ has 



no preference among them (wj wj' and Wj' wj). The first two cases are symmetrical 
and can be handled analogously. The last case follows from the first case because it has stronger 
assumptions. We prove the first case and assume that man prefers woman Wj to woman Wj> . 
The rules (|4]) imply the presence of a rule manpropose{mi,Wj') notaccept{mi,Wj), . . . 

and this is the only rule which can make manprcypose{mi,Wji) true (the only rule with this 
literal in the head). However, since accept{mi,Wj) is also in the answer set, this rule has a false 
body so manpropose{mi,Wj/) can never be in /. Consequently accept{mi,Wji) can never be in 
/ since the only rule with this literal in the head is of the form ([Ij and this body can never be 
true, which leads to a contradiction. 

Secondly assume that accept {mi, Wj) and accept{mi,mi) are both in an answer set / of V. 
Again step 1 implies that j G acceptable\,j . Because of the rules ^ V will contain the rule 
accept{mi,mi) ^ . . . ,notaccept{mi,Wj), . . .. An analogous reasoning as above implies that 
since accept{mi, Wj) is in the answer set /, accept^nii, mi) can never be in /. 

3. For every man mi, in every answer set I ofV exactly one of the following conditions is satisfied: 

(a) there exists a woman wj such that accept{mi,Wj) G /, 

(b) accept{mi,mi) e /, 

and similarly for every woman Wi . 

Suppose / is an arbitrary answer set of V and mi is an arbitrary man. We already know from 
step 2 that a man cannot be paired to a woman while being single, so both possibilities are 
disjoint. So suppose there is no woman Wj such that acceptijUi, Wj) is in /. V will contain the 
rule ([2]). Because of our assumptions and the definition of the reduct, this rule will be reduced 
to accept{mi,mi) ■<— , and so accept(mi,mi) will be in /. 

4. For an arbitrary answer set I of V the previous steps imply that / produces a set of marriages 
without blocking individuals. Weak stability also demands the absence of blocking pairs. Sup- 
pose by contradiction that there is a blocking pair (m^, Wj), implying that there exist i ^ i' and 
j ^ j' such that accept{mi,Wj') G I and accept{mii ,Wj) G / while Wj Wj' and mi <^ mi'. 
The rules of the form ([1]), the only ones with the literals accept(mi,Wji) and accept{mii , w j) 
in the head, imply that literals manpropose{mi,Wji) and womanpropose{mii ,Wj) should be 
in /. But since wj Wji and because of the form of the rules ([3]) there are fewer con- 
ditions to be fulfilled for manpropose{mi,Wj) to be in / than for manpropose{mi,Wjr) to 
be in /. So manpropose{mi,Wj) should be in / as well. A similar reasoning implies that 
womanpropose{mi, Wj) should be in /. But now the rules of the form ([l} imply that accept{mi, Wj) 
should be in /, contradicting step 2 since accept{mi,Wji) and accept{mii , w j) are already in /. 

□ 

Proposition 2 

Let {Sm, Sw) be an instance of the SMP with unacceptability and ties, and let V he the corresponding 
ASP program. If {{xi,yi), (xk,yk)} is a weakly stable matching for {Sm,Sw) then V has the 
following answer set I: 

{manpropose{xi, y) \ i e {1, . . . ,k},Xi e M, y y,)} 
U{womanpropose{x, yi)\i ^ {1, . . . , fc}, yi izW,x <^ Xi} 



U{accept{xi, yi) | i £ {1, . . . , k}} 
U{manpropose{xi, j/i) | « G {1, . . . , k}, Xi ^ y,} 
U{womanpropose{xi,yi) \ i € {1, . . . ,k},Xi ^ yi} 

Proof. Suppose we have a stable set of marriages S — {{xi,yi), . . . , {xk,yk)}, implying that every 
yi is an acceptable partner of Xi and the other way around. The rules of the form ([T]) do not alter 
when forming the reduct, but the other rules do as those contain naf-literals. Notice first that the 
stability of S implies that there cannot be an unmarried couple {m,w), with m a man and w a 
woman, such that manproposelm, w) is in / and womanpropose{m, w) is in I. By definition of / this 
would mean that they both strictly prefer each other to their current partner in S. This means they 
would form a blocking pair, but since S was stable, that is impossible. So the rules of the form ([T]) 
will be applied exactly for married couples (rui, Wj), since by definition of / manpropose{'mi, Wj) and 
womanpropose{mi,Wj) are both in / under these conditions. For other cases the rule will also be 
fulfilled since the body will be false. This reasoning implies that the unique minimal model of the 
reduct w.r.t. / should indeed contain accept{nii,Wj) for every married couple {mi,Wj) in S. Since S 
is a stable set of marriages, every person is either married or single. If a man rrii is single, there will 
be no other literal of the form accept{mi, .) in /, so rule ([2|) will reduce to a fact accept{mi, rrii) <— , 
which is obviously fulfilled by /. Similarly if a woman Wj is single. Any other rule of the form ([2|) or 
([3]) is deleted because rm or Wj is not single in that case and thus there is some literal of the form 
accept(mi,w) for some woman w and some literal of the form accept{m,Wj) for some man m in /, 
falsifying the body of the rules. If rrii is single, then accept{mi, rrii) is in / and this is the only literal 
of the form accept{mi, .) in J, so the rules of the form ^ will all be reduced to facts. The rule heads 
of these facts should be in the minimal model of the reduct and are indeed in I since the women w for 
which manpropose{mi, w) is in / are exactly those who are strictly preferred to staying single. The 
rules of the form Q for women Wj in neutral\j will all be deleted in this case, because accept{mi, rrii) 
is in /. If man mi is married to a certain woman Wj in the stable set S then the rules of the form Q 
will reduce to facts of the form manpropose{mi, w) ^ for every woman w who is strictly preferred to 
Wj and will be deleted for every other woman appearing in the head, because those rules will contain 
not accept(mi, wj) in the body. Again / contains these facts by definition, as the minimal model of the 
reduct should. We can use an analogous reasoning for the women. So the presence of the literals of 
the form manpropose{., .), womanpropose{., .) and accept{., .) in / is required in the unique minimal 
model of the reduct w.r.t. /. We have proved that every literal in I should be the minimal model of 
the reduct and that every rule of the reduct is fulfilled by /, implying that / is an answer set of 7-". □ 

In [18| it is shown that the decision problem 'is the pair (m, w) stable?' for a given SMP instance 
with unacceptablity and ties is an NP-complete problem, even in the absence of unacceptability. A 
pair {m,w) is stable if there exists a stable set that contains {m,w). It is straightforward to see that 
we can reformulate this decision problem as 'does there exist an answer set of the induced normal 
ASP program V which contains the literal accept{m^w)T (i.e. brave reasoning), which is known to 
be an NP-complete problem [I] . So our model forms a suitable framework for these kind of decision 
problems concerning the SMP. 



4 Selecting Preferred Stable Sets 



4.1 Notions of Optimality of Stable Sets 

When several stable matcliings can be found for an instance of the SMP, some may be more interesting 
than others. The stable set found by the G-S algorithm is M-optimal [35], i.e. every man likes this set 
at least as well as any other stable set. Exchanging the roles of men and women in the G-S algorithm 
yields a W- optimal stable set [T, optimal from the point of view of the women. 

While some applications may require us to favour either the men or the women, in others it makes 
more sense to treat both parties equally. To formalize some commonly considered notions of fairness 
and optimality w.r.t. the SMP, we define the cost Cx{S) of a stable set S to an individual x, where 
Cx{S) — k a X has been matched with his or her k*^ preferred partner. More precisely, for x — rrii a 
man, we define c„n{S) — \{z : z ?/}| + 1 where y is the partner of x in S] for x = Wj a woman, Cx 
is defined analogously. So in case of ties we assign the same list position to equally preferred partners, 
as illustrated in Example [S] 

Example 5 

Let X — mi be a man with preference list a]^j = ({1}, {2, 3}, {4}) then wi as partner of x in some set 
of marriages S would yield Cx{S) ~ 1, W2 and W3 yield Cx{S) = 2 and W4 yields Cx{S) = 4. If mi 
would be single in S, then the cost Cx{S) is 4, since mi prefers women wi,W2 and ^3 to being single, 
but is indifferent between being paired to W4 or staying single. 

Definition 5 

For S a set of marriages, 

• the sex-equalness cost is defined as Csexeq{S) = \ J2xeM {S)\, 

• the egalitarian cost is defined as Cweight{S) — J2xeMiiw '^^(^^f 

• the regret cost is defined as Cregret{S) = UY&yLxszMvjw Cx{S) , and 

• the cardinality cost is defined as Csingies{S) = \{z : {z,z) G S}\. 

S is a sex-equal stable set iff S is a stable set with minimal sex-equalness cost. Similarly, S is an 
egalitarian (resp. minmum regret, maximum cardinality) stable set iff S is a stable set with minimal 
egalitarian (resp. regret or cardinality) cost. 

A sex- equal stable set assigns an equal importance to the preferences of the men and women. An 
egalitarian stable set is a stable set in which the preferences of every individual are considered to 
be equally important. In the use of an egalitarian stable set is proposed to optimally match 
virtual machines (VM) to servers in order to improve cloud computing by equalizing the importance 
of migration overhead in the data center network and VM migration performance. A minimum regret 
stable set is optimal for the person who is worst off. A maximal or minimal cardinality stable set 
is a stable set with resp. as few or as many singles as possible. Examples of practical applications 
include an efficient kidney exchange program [3T] and the National Resident Matching Progrand [l5] . 
Maximizing cardinality garantuees that as many donors as possible will get a compatible donor and 
as many medical graduates as possible will get a position. 



'www.nrmp.org 



Table [T] presents an overview of known complexity resultf0 concerning finding an optimal stable 
set. Typically the presence of ties leads to an increase of complexity. Manlove et al. |T71 [TS] proved 
that the problem of finding a maximum (or minimum) cardinality stable set for a given instance of the 
SMP with unacceptability and ties is NP-hard. Using this result, the problem of finding an egalitarian 
or minimum regret stable matching for a given SMP instance with ties is proved to be NP-hard [18] . 
even if the ties occur on one side only and each tie is of length 2 (i.e. each set in a preference list has 
size at most 2). If there are no ties, the problem of finding an egalitarian or minimum regret stable set 
is solvable in polynomial time [131 [TO] • Since all stable sets consist of n couples in the classical SMP, 
the G-S algorithm trivially finds a maximum (or minimum) cardinality [7]. For an SMP instance 
with unacceptability the number of couples in a stable set is constant [8], so finding a maximum 
cardinality stable set reduces to finding a stable set, which is known to be solvable in polynomial 
time. Surprisingly, finding a sex-equal stable set for a classical SMP instance is NP-hard [16], even if 
the preference lists are bound in length by 3 [20] • 



Table 1: Literature complexity results for finding an optimal stable set 
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Between brackets we mention in Table [T] the complexity of an algorithm that finds an optimal 
stable set if one exists, in function of the number of men n. To the best of our knowledge, the 
only exact algorithm tackling an NP-hard problem from Table [1] finds a sex-equal stable set for an 
SMP instance in which the strict preference lists of men and/or women are bounded in length by 
a constant [20]. To the best of our knowledge, no exact implementations exist to find an optimal 
stable set for an SMP instance with ties, regardless of the presence of unacceptability and regardless 
which notion of optimality from Table [1] is used. Our approach yields an exact implementation of all 
problems mentioned in Table [TJ 

4.2 Finding Optimal Stable Sets using Disjunctive ASP 

As we discuss next, we can extend our ASP encoding of the SMP such that the optimal stable sets 
correspond to the answer sets of an associated ASP program. In particular, we use a saturation 
technique [5] [T] to filter non-optimal answer sets. Intuitively, the idea is to create a program with 
3 components: (i) a first part describing the solution candidates, (ii) a second part also describing 
the solution candidates since comparison of solutions requires multiple solution candidates within the 
same answer set whereas the first part in itself produces one solution per answer set, (iii) a third 
part comparing the solutions described in the first two parts and selecting the preferred solutions by 
saturation. It is known that the presence of negation- as-f allure can cause problems when applying 
saturation. Therefore, we use a SAT encoding [T5] of the ASP program in Definition U and define a 



■^Throughout this paper we assume that P 7^ NP. 



disjunctive naf-free ASP program in Definition |5] which selects particular models of the SAT problem. 
Notice that our original normal program is absolutely tight, i.e. there is no finite sequence li,l2, ■ ■ ■ 
of literals such that for every i there is a program rule for which li^i is a positive body literal and 
li is in the head ^ . We use the completion and a translation of our ASP program to SAT to derive 
Definition [S] The completion of a normal ASP program is a set of prepositional formulas. For every 
atom a with a bodyi (i G {1, . . . , fc}) all the program rules with head a, the propositional formula 
a = body'i V ... V body'/^ is in the completion of that program. If an atoms a of the program docs not 
occur in any rule head, than a = _L is in the completion of the program. Similarly the completion of the 
program contains the propositional formula _L = body[ V ... V body'^ with with ^ bodyi [i G {!,...,/}) 
all the program constraints. For every z, body'^ is the conjunction of literals derived from bodyi by 
replacing every occurence of 'noi' with Because our program is absolutely tight, we know that 
the completion will correspond to it [6]. Applied to the induced normal ASP program in Definition 
I3J the completion becomes: 

{acceptlrrii, Wj) = manproposeirrii, Wj) A womanpropose{mi, Wj) \ i G {1, . . . , n},j G {1, . . . ,p}} 
U{accept{mi,mi) = ^ -iaccept{mi,Wk) \i ^ {1, ■ ■ ■ ,n}} 

kGacceptable\j 

U{accept{wj , Wj) = ^ ^accept{mk, Wj) | j G {1, . . . ,p}} 

k^acceptable^^ 

\j{manpropose[mi, Wj) = ^ ^accept{mi, x) | i G {1, . . . , n}^j G acceptable\,f} 
U{womanpropose{mi, Wj) = -^accept{x, Wj) \ j e {1, . . . ,p}, i G acceptableiy} 

x<^Y ,x^mi 

U{manpropose{'mi,Wj) = ± |i G {!,..., ri},j G unacceptable\.j} 

U{womanpropose{mi, Wj) = 1. \jE {1, . . . i G unacceptable-^^} 

Using the formulas of the completion corresponding to the normal ASP program in Definition |4l we 
can define a corresponding disjunctive ASP program without negation-as-failure. Lemma |3] follows 
form the fact that the completion corresponds to the original program [5]. 

Definition 6 (Induced disj. naf-free ASP program) 

The disjunctive naf-free ASP program Vdisj induced by an SMP instance {Sm, Sw) with unacceptability 
and ties contains the following rules for i G {1, . . . ,ri}, j G {1, . . . ,p} : 

-'accept{mi, Wj) V manpropose{mi, Wj) <— 
-iaccept{mi, Wj) V womanpropose{mi, Wj) <— 
accept{mi,Wj) V -'manpropose{mi,Wj) V -iwomanpropose{mi,Wj) 
For every i G {1, . . . , n} , I G unacceptable\j , j G acceptable\j, x <™i Wj,x ^ Wj Vdisj contains: 

\J accept{mi, Wk) V accept{mi, mi) ^ 

kGacceptable^j^j 

-^accept{mi, rui) V -iaccept{mi, Wj) 
-imanproposelmijWj) V -'accept(mi, x) 



\J accept{mi, x) V manpropose[mi, wj) ■<— 

-^manpropose{mi,wi) Br- 
unei symmetrical for j £ {1, . . . ,p} and womanpropose. 

Lemma 3 

Let V he the normal ASP program from Definition^ and Vdisj the disjunctive ASP program from 
Definition\^ It holds that for any answer set lofV there exists an answer set Idisj of Vdisj such that 
the atoms of I and Idisj coincide. Conversely for any answer set Idisj of Vdisj there exists an answer 
set I of V such that the atoms of I and Idisj coincide. 



4.3 ASP Program to Select Optimal Solutions 

Let {Sm,Sw) be an SMP instance with unacceptability and ties, with Sm ~ Wmt--t'^m} ^-nd 
Sw = ■ ■ ■ j^w}' ^^'^ 'Pnorm be the induced normal ASP program from Definition |31 Our 

technique for extending this program to a program that can respectively optimize for the sex-equalness, 
egalitarian, minimum regret and maximum cardinality criterion is in each case very similar. We start 
by explaining it for the case of sex-equalness. Our first step is to add a set of rules that compute the 
sex-equalness cost of a set of marriages. For every man and every woman Wj such that j G (7\j(k) 



we use the following rule to determine the cost for if Wj would be his partner: 

mancost{i, k) accept{mi,Wj) (6) 
and similarly for every Wj and every rrii such that i £ CT^(fc): 

womancost{j,k) <^ accept{mi,Wj) (7) 
We also use the following rules with i ranging from 1 to n and j from 1 to p: 

mancost{i, |cr\,/|) <— accept{mi,mi) (8) 
womancost{j, \cr^,\) <— accept{wj,Wj) (9) 
manweight(Z) ^ f/=sum{B, A : mancost{A, B)} = Z, jj^int{Z) (10) 
womanweight{Z) ^ ff:sum{B, A : womancost{A, B)} = Z, ^int{Z) (11) 
sexeq{Z) ^ manweight{X), womanweight{Y) , Z ^ X — Y 

sexeq{Z) <— manweight{X), womanweight{Y) , Z = Y — X (12) 



Rules dS]) and © state staying single leads to the highest cost. Rule pUj) determines the sum of 
the male costs^ and similarly pT|) determines the sum of the female costs. According to Definition 
[S]the absolute difference of these values yields the sex-equalness cost, as determined by rules p^ . 
Since numeric variables are restricted to positive integers in DLV, we omit conditions as 'X > Y^ or 
'X < Y\ The program Vnorm extended with rules (jH) - (fT2|) is denoted Vl^^'^'^ . We construct a 
program Vsexeq, composed by subprograms, that selects optimal solutions. Let V'^i^j be the disjunctive 
naf-free ASP program, induced by the same SMP instance, in which a prime symbol is added to all 
literal names (e.g. accept becomes accept'). Define a new program V'J^^^^'^'^ with all the rules of P^j^, in 



^i^sum, ^max, #int and ^count are DLV aggregate functions. The 'A' mentioned as variable in #sum indicates 
that a cost must be included for every person (otherwise the cost is included only once when persons have the same 
cost). 



which every occurrence of -latom is changed into natom for every atom atom, i.e. replace all negation 
symbols by a prefix 'n'. For every occurring atom atom in V'^^i'^'^, add the following rule to exclude 
non-consistent solutiont|f|: 

sat <— atom, natom (13) 

Finally add rules ([6]) - (fT2|) with prime symbols to the literal names to V'fJt^'^'^ but replace rule ([TOl) 
and rule ([TT]) by: 

mansum{n, X) <— mancost{n, X) 

mansum{J, Z) ^ mansum{I , X), mancost{J, Y), Z = X + Y, jfsucc{J, I) 
manweight{Z) ^ mansum{l, Z) 
'womansum{p, X) ^ womancost{p, X) 

womansum{J, Z) ^ womansum{I , X), womancost{J, Y), Z = X + Y, ^succ{J, I) 

womanweight{Z) ^ womansum{l, Z) (14) 

The DLV aggregate function ^succ{J, I) is true whenever J + 1 = /. The reason we replace the rules 
with the aggregate function ^max by these rules is to make sure the saturation happens correct. When 
saturation is used, the DLV aggregate function #max, jj^sum and ^count would not yield the right 
criteriumvalues. Moreover, DLV does not accept these aggregate function in saturation because of the 
cyclic dependency of literals within the aggragate functions created by the rules for saturation. These 
adjusted rules, however, will do the job because of the successive way they compute the criteriumvalues. 
This becomes more clear in the proof of Proposition^ We define the ASP program Vsexeq as the union 
of Vext'^'^ J 'PexT'^'^ and Vsat- The ASP program Vsat contains the following rules to select minimal 
solutions based on sex-equalness: 

sat sexeq{X),sexeq'{Y), X <Y (15) 
■f- not sat (16) 
mancost'{X,Y) sat,manargcost'i{X),manargcost'2{Y) 
womancost' {X,Y) ^ sat,womanargcost'i{X),womanargcost'2{Y) (17) 
manpropose (X,Y) sat,man{X),'woman(Y) 
womanpropose' {X,Y) <— sat,man[X),woman{Y) 
accept' {X, X) <— sat,man{X) 
accept' {X, X) sat,'woman{X) 

accept' {X,Y) <— sat,man{X),woman{Y) (18) 

and analogous to p8)) a set of rules with prefix 'n' for the head predicates. Finally we add the fact^ 
manargcost'i{l..n) manargcost'2{l-.ip + 1)) womanargcost'i(\..p) womanargcost'2{l..{n + 
1)) <— , man(x) •<— for every man x and woman{x) •(— for every woman x to Vsat- Intuitively the 
rules of Vsat express the key idea of saturation. First every answer set is forced to contain the atom 
sat by rule ([TB)) . Then the rules ((T7]) - and the facts make sure that any answer set should 
contain all possible literals with a prime symbol that occur in Vsexeq- Rule (jl5l) will establish that 
only optimal solutions will correspond to minimal models and thus lead to answer sets. For any 

^For instance, sat <— accept' (mi, wi),naccept' {m\,wi) 

^The rule manargcost'^{l..ri) •<— is DLV-syntcLX for the n facts manargcost'^{l) <— , . . . ,manargcost'-^{n) <— . 



non-optimal solution, the corresponding interpretation containing sat will never be a minimal model 
of the reduct. It is formally proved in Proposition 2] below that Vsexeq produces exactly the stable 
matchings with minimal sex-equalness cost. 

Furthermore, only small adjustments to Vsexeq are needed to create programs Vweight, Vregret, 
and Vsingies that resp. produce egalitarian, minimum regret and maximum cardinality stable sets. 
Indeed, the ASP program Vweight can easily be defined as Vsexeq in which the predicates sexeq and 
sexeq' are resp. replaced by weight and weight' and the rules ([T^ are replaced by ([T5|. determining 
the egalitarian cost of Definition [S] as the sum of the male and female costs: 

weight{Z) ^ manweight{X), womanweight{Y) , Z — X + Y (19) 

Similarly the ASP program Vregret is defined as Vsexeq in which the predicates sexeq and sexeq' 
are resp. replaced by regret and regret' and rules (1101) - (jl2p are replaced by the following rules: 

manregret{Z) ^ jj=max{B : mancost{A, B)} = Z, ^int{Z) (20) 

womanregret(Z) <— =ffmax{B : womancost{A, B)} = Z,^int{Z) (21) 

regret{X) ^ manregret{X),womanregret{Y), X > Y 

regret{Y) ^ manregret{X),womanregret{Y), X < Y (22) 

Rule (|20|) determines the regret cost but only for the men. Similarly (|21|) determines the regret cost for 
the women. The regret cost as defined in Definition [5] is the maximum of these two values, determined 
by the rules in (f^ . Again we adjust rules (PO]) and (PTI) for the program part V'J^f^"^* by replacing 
them with a successively computing variant: 

manmax{n, X) -f- mancost{n, X) 

manmax{J, X) 4— manmax{I , A"), mancost{ J, Y), X >Y, ^succ{J, I) 
manmax{ J, Y) <— manmax{I , A), mancost{ J, Y), A < Y, ^succ{J, I) 
manregret{Z) ^ manmax{l, Z) 
womanmax{p, A) <— womancost(p, A) 

womanmax{ J, A) womanmax{I , A), womancost{ J, Y), A > Y, ^succ{J, I) 
womanmax{J, Y) <— womanmax{I , A), womancost{ J, Y), A < Y, jfsucc{J, I) 
womanregret{Z) -f- womanmax{l, Z) (23) 

Finally we define the ASP program Vsingies as Vsexeq in which the predicates sexeq and sexeq' are 
resp. replaced by singles and singles'. Furthermore we replace rules ([6|) - (|12|) by (|24|) . determining 
the number of singles: 

singles{Z) #count{B : accept{B, B)} = Z, #int{Z) (24) 

This time we adjust rule ([24|) for the program part V'^^f^^^'' as follows: 
single{p + i, 1) accept{mi, rrii) 
single(p + i, 0) <— naccept{mi, rrii) 
single{j, 1) ^ accept{wj,Wj) 
single{j,0) ^ naccept{wj,Wj) 
singlesum{n + p. A) single{n + p, A) 

singlesum{J, Z) •<— singlesum{I , A), single{J, Y), Z = X + Y, =fl=succ{ J, I) 



singles{Z) ^ sat, singlesum{l, Z) (25) 

Example 6 

We reconsider Example |31 This SMP instance had 3 stable sets of marriages: 

• Si = {accept{mi,'W3),accept{m2,'Wi),accept{w2,W2)}, 

• S2 — {accept{mi,W2),accept{m2,'Wi),accept{w3,W3)}, 

• 83 = {accept{rni,'Wi),accept{rn2,m2),accept{w2,'W2),accept{w3,W3)}. 

It is easy to compute the respective regret costs as Cregret{Si) = 2 and Cregret{S2) = CregretiSs) — 3. 
The corresponding program selecting this minimum regret stable set is the program consisting of the 
rules in Example 2] in addition with: 
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mancost' {1, 2 
mancost' {1, 2 
mancost' {1, 3 
m.ancost' {2, 2 
mancost' {2, 1 
mancost' {2, 2 
womancost'{l, 1 
womancost' {1, 1 
womancost' (1, 2 
womancost' {2, 1 
womancost' {2, 2 
womancost' {3, 2 
womancost' {3, 1 
womancost' {3, 3 
manmax' {2, X 
manmax' {J, X 
manmax' {J, X 
manregret' {Z 
womanmax' {2, X 
womanmax' {J, X 
womanmax' (J, X 
womanregret' {Z 
regret' {X 
regret' {Y 
sat 

manargcost'i{1..2 
manargcost'2 ( 1 . . 4 
womanargcost'i ( 1 . . 3 
womanargcost'i (1..3 
mancost' {X,Y 
womancost' {X , Y 
manpropose' {X , Y 
nmanpropose' {X, Y 



accept' {X, X), naccept'{X, X), woman{X) 

accept' {mi, wi) 

accept' {mi, W2) 

accept' {mi, wz) 

accept' {mi, mi) 

accept {m,2, wi) 

accept' (m2, W2) 

accept' {m2, m2) 

accept' {mi, wi) 

accept' {m2,wi) 

accept' {wi, wi) 

accept' {mi, W2) 

accept' {w 2, w 2) 

accept' {mi, ws) 

accept' {m2, w^) 

accept' {w 3, w 3) 

mancost' {2, X) 

manmax' {I, X), mancost' {J, Y),X Y, ^succ{J, I) 
manmax' {I, X), mancost' {J, Y),X >= Y, =ffsucc{J, I) 
manmax' {1, Z) 
womancost' {2, X) 

womanmax' {I, X) , womancost' {J, Y),X >= Y, ^succ{J, I) 
womanmMx {I , X), womancost' {J, Y), X >= Y, ^succ{J, I) 

womanmax' {1, Z) 

manregret' {X), womanregret' {Y),X > Y 
manregret' {X),womanregret'{Y),X <= Y 
regret{X), regret' {Y),X <= Y 
not sat 



i— sat, manargcost'i {X), manargcost'2 0^) 

<— sat, womanargcost'i{X), womanargcost'2{Y) 

sat, man{X), womaniY) 
•(— sat, man{X), woman{Y) 



womanpropose' {X,Y) ^ sat,man{X),woman{Y) 

nwomanpropose {X,Y) sat,man{X),woman{Y) 

accept' {X, Y) sat, man{X), woman{Y) 

accept' {X, X) <r- sat,man{X) 

accept' (X,X) -s— sat,woman{X) 

naccept' {X ,Y) ^ sat,man{X),woman{Y) 

naccept' {X, X) <r- sat,man{X) 

naccept' {X, X) ^ sat,woman{X) 

Letting DLV compute the unique answer set of this disjunctive ASP program and filtering it to 
the hterals accept and regret, yields {accept{m2,wi), accept (m 1,1113), accept{w2,W2), regret(2)}, 
corresponding exactly to the minimum regret stable set of the SMP instance and the corresponding 
regret cost. 

Proposition 4 

Let the criterion crit be an element of {sexeq, weight, regret, singles}. For every answer set I of 
the program Vcrit induced by an SMP instance with unacceptability and ties the set Si — {{m,w) 
accept(m, w) G /} forms an optimal stable set of marriages w.r.t. criterion crit and the optimal 
criterion value is given by the unique value vj for which critivj) £ I . Conversely for every optimal 
stable set S — {{xi,yi), .. ., {xk,yk)} with optimal criterion value v there exists an answer set I of 
Vcrit such that {{x,y) \ accept{x,y) G /} = {(xi,yi) \ i £ {l,...,fc}} and v is the unique value for 
which crit{v) £ /. 

Proof. Let {Sm, Sw) is an instance of the SMP with unacceptability and ties. 



Answer set => Optimal stable set Let / be an arbitrary answer set of Vcrit and let Si be as in 
the proposition. It is clear that the only rules in Vcrit that influence the literals of the form 
manpropose{., .), womanpropose{., .) and accept{., .) are the rules in Vnorm- Hence any answer set / 
of Vcrit should contain an answer set Inorm of Vnorm as a subset. Proposition [T] implies that Inorm 
corresponds to a stable set Si — {{m,w) \ accept{m,w) £ Inorm}- Moreover, the only literals of the 
form manpropose(., .), womanpropose{., .) and accept{., .) in / are those in Inorm, so 5*/ = {{m, w) 
accept{m,w) £ /}. If crit = sexeq, it is straightforward to see that the literals of the form accept{., .) 
in Inorm Uniquely determine which literals of the form mancost{., .), womancost{., .), manweight{.) , 
womanweight{.) and sexeq{.) should be in the answer set /. These literals do not occur in rules of 
Vcrit besides those in V^xt'^'^ ■ Notice that the rules which do contain these literals will imply that 
there will be just one literal of the form sexeq{.) in /, namely sexeq{v) with v the sex-equalness cost 
of Si. Analogous results can be derived for crit £ {weight, regret, singles} . It remains to be shown 
that Si is an optimal stable set. Suppose by contradiction that Si is not optimal, so there exists a 
stable set S* such that vi > v* , with v* the criterion value of S* to be minimized. We will prove that 
this implies that / cannot be an answer set of Vcrit , contradicting our initial assumption. 
Proposition [2] and Lemma [3] imply that there exists an interpretation I%igj of the ASP program Vdisj 
induced by Sw) that corresponds to the stable set S* . Moreover this interpretation is consistent, 
i.e. it will not contain atom and -latom for some atom atom. This implies that the interpretation 
^'disj defined as I^isj in which -latom is replaced by natom for every atom atom will falsify the body 
of the rules of the form ([T5)) of V'^^^t* . An analogous reasoning as above yields that the literals of the 



form accept' {.^ .) in I'^i^j uniquely determine which literals of the form mancost' {., .), womancost' .), 
mansum' {., .), womansum'l., .), manweight'{.), womanweight' {.) and sexeq'{.) should be in /^j^j. 
With those extra literals added to I'^^sj, I'^i^j satisfies all the rules of Vj^f*. Moreover, crit{v*) is the 
unique literal of the form crit{.) in I'^^^^y Notice that /^j^^- does not contain the atom sat. 
Define the interpretation J — Inorm U I'disj ■ From the previous argument it follows that J will satisfy 
every rule of Vll.f U V^^^^ since the predicates occurring in both programs do not overlap. Moreover 
J contains critivi) and crit'(v*) and these are the only literals of the form crit{.) or crit'{.). Since 
vj > V* the rules of the form ^T5\i will be satisfied by J since their body is always false. Call J' the 
set J U {a\{a ^) G Vsat}- Since J' does not contain sat, the rules of Vsat will all be satisfied by J', 
with exception of the rule ^ not sat. 

The rule of the form implies that / as answer set of Vcrit should contain sat. Now the set 
of rules ([T7)) - imply that / should also contain the literals mancost' {.,.), womancost' {., .) and 
manpropose' .), womanpropose' .), accept' .) with the corresponding literals prefixed by n for 
every possible argument stated by the facts in Vsat- The successively computing rules (ITU) resp. ([^ 
and ^ in V'^Zt \ by which we replaced rules ^ - ^ resp. ^ - ^ and (HH), garantuee that 
for every possible set of marriages and its corresponding criterium value c / will contain crit{c) and 
all associated intermediate results. E.g. for crit — sexeq, the rules will garantuee that / also contains 
mansum' {., .), manweighti.) , womansum{., .) and womanweight{.) for every argument that could 
occur in a model of "Pez:"*- □ This implies that I'^^^^j C /. We already reasoned in the beginning of the 
proof that Inorm Q I holds so it follows that J C J. Since the literals oi J' \ J are stated as facts of 
V^"/, they should be in /, hence J' C /. Moreover J' C I since sat £ / \ J'. 

We use the notation red{V, I) to denote the reduct of an ASP program P w.r.t. an interpretation /. 
There is no rule in Pex7* with negation-as-failure in the body, hence red{Vj^Jf*, I) = red{V'^^i^ , J') ~ 
'P'eZt*- We already reasoned that J' satifies all the rules of the latter. We also reasoned that / does not 
contain any other literals of the form accept (., .) than those who are also in Inorm, and by construction 
the same holds for J'. Hence red{V^"i,I) = redCP^J^f , J') and by construction J' satisfies all the 
rules of this reduct. It is clear that red{'Psat, I) is Vsat without the rule •(— not sat, since sat G /. 
Again we already argued that J' satisfies redCPsat, I)- Hence J' satisfies all the rules of red{Vcrit, I), 
implying that /, which strictly contains J', cannot be an answer set of Vcrit since it is not a minimal 
model of the negation- free ASP program red{Vcrit, I) [S]- 



Optimal stable set => Answer set Let {Sm,Sw) be an instance of the SMP with unacceptability 
and ties and let S = {{xi, yi), . . . , (xfe, yk)} be an optimal stable set with optimal criterion value v. 
To see that the second part of the proposition holds it suffices to verify that the following interpre- 
tation / is an answer set of Vcrit, with the notation Pxi{y) as the index a for which y G (j\,j{a) if 
Xi = mi and symmetrically Py. (x) as the index a for which x G cr^^ (a) if yi = wi' . If Xi — yi we set 
P'Xiiyi) = Pyi{xi) = WXil if Xi is a man and \a\y \ otherwise. So let / be given by: 

/ = /i U /2 



^Notice that this would not be the case if we use the original rules with ^sum, ^max and ^count in P^^"', since 
these rules would lead to only one value cm for which e.g. manweight{cM) should be in I, and similarly only one value 
cw for which womanweight{cw) should be in /. Consequently there would be only one value c such that crit(c) should 
be in /. This value would not necessarily correspond to v* and so we would not be able to conclude that I'^^^j C /. 
Moreover DLV does not allow the use of these rules because of the cyclic dependency of literals they would create, 
involving the variables in the aggregate functions. 



with 

Ii ={accept{xi, yi)\i €z {1, . . . , k}}{crit{v)} U {sat} 

U{womanpropose{xi,yi) \ Xi ^ yi}{manpropose{xi,yi)\xi 7^ j/j} 
U{manpropose{xi,y)\i e {l,...,k},Xi =mi,3a< Pxiivd'-V ^ <^m{(^)} 
U{womanpropose{x,yi)\i e {l,...,k},yi = wi',3a < Py^{xi):x G cr^(a)} 
\j{mancost{l, Pxi{yi)) \ crit ^ singles, « G {1, . . . , k}, Xi = mi} 
\j{'womancost(Py. (x^), I') \ crit ^ singles, i G {1, . . . , k}, yi = wi'} 
U{manweight{cM{S)) \crit G {sexeq, weight}} 
U{womanweight{cwiS)) | crit G {sexeq, weight}} 
U{manregret{cregret,M{S)) \ crit = regret} 
U{womanregret{cregret,w{S)) \ crit = regret} 

and 

I2 ={manargcost'i{z) \z G {l,...,n}} 
(j{rnanargcost'2{z) \ z G {1, . . . ,p + 1}} 
yj {womanarg CO st'i{z) \z G {l,...,p}} 
U{womanargcost'2{z) \ z G {1, . . . ,n + 1}} 

U{man{x) \ x G M} U {woman{x) \ x G W} (26) 
U{mancost' | crit / singles, i G {1, . . . , n},j G {1, . . . ,p + 1}} 

U{womancost' (j , i) \ crit 7^ singles, i G {1, . . . , n + 1}, i G {1, . . . ,p}} (27) 
U{manpropose' {x, y) | x G M, y G W} U {womanproposeJ {x, y) | x G M, ?/ G W^} 
U{accept'(a;, y) | x G M, y G VF} U {accept' {x, x) | x G M U W^} 
U{nmanpropose'{x,y) \ x G M,y e. W} 
U{nwomanpropose'{x,y) \ x G M,y G W} 

Ll{naccept'{x, y) | x G M, y G VK} U {naccept'{x, x)\x G MUW} (28) 
U{crit' (val) \ val G arg(crit)} 

U{single'{i, j) \ crit = singles, i G {1, . . . ,n + p}, j G {0, 1}} 

U{singlesum'{i, j) | crit — singles, i G {1, . . . ,n + p},j G {1, . . . , {n + p — i + 1)}} 
U{mansum'{i,j) \ crit G {sexeq, weight}, i G {1,. . ., n},j G {n — i + 1, . . . , {n — i + l)(p + 1)}} 
U{womansum'{j, i) \ crit G {sexeq, weight}, j G {1, . . . ,p},i G {p — j + 1, . . . , {p — i + l){n + 1)}} 
U{manweight' (z) \ crit e {sexeq, weight}, z G {n, . . . , n{p + 1)}} 
U{wcymanweight' (z) \ crit G {sexeq, weight}, z G {p, ■ ■ ■ ,p(n + 1)}} 
yj {manmax {i, j) \ crit = regret, i G {1, . . . , n},j G {1, . . . ,p + 1}} 
U{womanmax'{j, i) \ crit = regret, j G {1, . . . ,p}, i G {1, . . . , n + 1} 
U{manregret' (z) \ crit = regret, z G {1, . . . ,p+ 1}} 

U{womanregret' {z) \ crit = regret, z G {1, . . . ,n + 1}} (29) 



The notation arg(c) stands for the possible values the criterion can take within this problem instance: 

• crit = sexeq arg(crit) = {0, . . . , max(np + n — p, np + p — n)} , 

• crit = weight arg(crit) = {n + p, . . . , 2np + p + n}, 

• crit = regret arg(crit) = {1, . . . , max(p, n) + 1)}, 

• crit — singles arg(cri<) = {0, . . . ,ji + p}. 

To verify wether this interpretation is an answer set of Vcrit, we should compute the reduct w.r.t. / 
and check wether J is a minimal model of the reduct [5] . It can readily be checked that / satisfies all 
the rules of red{Pcrit, I)- It remains te be shown that there is no strict subset of / with satisfies all 
the rules. First of all all the facts of Vcrit must be in the minimal model of the reduct, explaining 
why the sets of literals should be in /. The only rules with negation-as- failure are part of V^xf- 
As in the previous part of the proof, it is straightforward to see that Ii is the unique minimal model 
of the reduct of "Pg^f w.r.t. /, considering that the literals in I2 don't occur in V^"/. So any minimal 
model of red{Vcrit , I) must contain Ii . 

The key rule which makes sure that / is a minimal model of the reduct is (jisp . The rules (fT3|) imply 
that for each model of red{Vcrit, I) that does not contain sat, the literals of V'fjf* in that model 
will correspond to a stable set of the SMP instance. In that case rule (|15l) will have a true body, 
since S is optimal, implying that sat should have been in the model. And the presence of sat in any 
minimal model implies the presence of the set of literals (I29p in any minimal model of the reduct. 
This can be seen with the following reasoning. Due to the presence of the facts ([^5]) and sat in 
any minimal model of the reduct, rules (1171) imply the presence of the literals (|27p in any minimal 
model. For the same reason rules ([T8| imply that the literals ((28)) should be in any minimal model 
of red{'Pcrit, I)- For crit = sexeq the presence of the literals of the form (j27p in any minimal model 
of the reduct together with rules (fM)) imply that mansum'{i,j) should be in any minimal model for 
every i € {1, . . . ,n} and j G {n — i + 1, . . . , {n — i + l){p + 1)}: for « = n the first rule of ([M]) implies 
that mansum'(n,x) is in any minimal model for every x such that manargcost'2{x) is in it, i.e. any 
X G {1, . . . Now the second rule of (1141) implies that mansum! {n— 1, x) is in any minimal model 

for every x-\-y such that manargcost'2{x) and mansum' (n, y) are in it, i.e. any x-\-y G {2, . . . , 2(p+l)}. 
If we continue like this, it is straightforward that every literal of the form mansum' .) of I2 should 
be in any minimal model. The third rule of (|14[) now implies that manweight' {x) should be in 
any minimal model for every x such that mansum' {l,x) is in it, i.e. x G {n, . . . , n(p + 1)}. The 
same reasoning can be repeated for the literals womansum' and womanweight' . At this point rules 
(jl2[) imply that sexeq' {\x — y\) should be in any minimal model which contains manweight' {x) and 
womanweight' {y) . Notice that only one of the two rules in (|12[) will apply for every x and y since 
the numerical variables in DLV are positive. Considering the arguments for which manweight' and 
womanweight' should be in any minimal model, it follows that sexeq' (x) should be in any minimal 
model for every x G {0, . . . , max(p(n + 1) — n, n{p+ 1) — p)}, which is exactly arg(crii). For the other 
criteria, an analogous reasoning shows that the presence of all literals of I2 is required in any minimal 
model of the reduct. 

Considering the fact that we have proved that all literals of / should be in any minimal model of the 
reduct and / fulfils all the rules of the reduct, we know that / is a minimal model of the reduct and 
thus an answer set of Vcrit- 

□ 



If we delete from Vsexeq the rules ([TT|) - and replace rule (IT5|) by the rule sat ^ ■manweight{X) , 
manweight'{Y),X < Y, then we obtain the M-optimal stable sets. Analogously we can obtain the 
W-optimal stable sets. 

If a criterion is to be maximized, the symbol < in rule (jlSp is simply replaced by >. E.g. for 
crit = singles we will get minimum cardinality stable sets. 

5 Complexity and Future Work 

The NP-complete decision problem 'does there exist a stable set with cardinality > k (resp. < k) for 
an SMP instance with unacceptability and ties with k a positive integer?' |17[|18j has practical impor- 
tance, e.g. in the National Resident Matching Program [TS|. If we add a rule sat ^ singles{X), X < 
{n+p — 2k) to the extended induced program V^^^''^" defined in Subsect. [4.3[ then this problem can be 
formulated as 'does there exist an answer set of the normal ASP program V^J^^'''^'^ which contains the 
literal sail' (i.e. brave reasoning), another NP-complete problem [T|. So our model forms a suitable 
framework for these kind of decision problems concerning optiniality of stable sets in the SMP. 

Notice that the complexity of this kind of decision problem and the one mentioned in the last 
paragraph of Subsect. |3] are a good indication how hard it is to find an (optimal) stable set, as 
opposed to the problems 'does there exist an (optimal) stable set?', which tell us how hard it is to 
know whether there exists a solution but not necessarily how hard it is to find one. 

Combining these problems leads to a new decision problem: 'is the pair (m, w) optimally stable 
for an instance of the SMP with unacceptability and ties?'. We define an optimally stable pair as a 
pair (m, w) for which there exists an optimal stable set in which m and w are matched. As far as 
we know this problem has not been studied yet, although it could be useful in practice, for instance 
if one wants to find a maximum cardinality matching but also wants to prioritize some couple or a 
person. Optimality is still desirable, because it ensures the others from not being put too much at a 
disadvantage. For instance in the kidney exchange problem, in which kidney patients with a willing 
but incompatible donor try to interchange each other's donors to get a transplant, this is a realistic 
situation: if two patients with intercompatible donors urgently need a transplant, they should get 
priority, but of course we still want to match as many patients to donors as possible. Considering the 
complexity of the separate decision problems, the combined problem might have a higher complexity, 
perhaps corresponding to the -complexity of our grounded disjunctive normal ASP program with 
aggregate functions [H E] . It should be noticed however that the addition of constraints not necessarily 
increases complexity and a precise classification of complexity is desirable. 

6 Conclusion 

We formalized and solved different variants of the SMP using ASP programs, which can easily be 
adapted to yet other variants. Moreover we applied saturation to compute optimal stable sets, with 
the advantage that these programs can be handled with the efficient off-the-shelf ASP solver DLV. 
To the best of our knowledge, our encoding offers the first exact implementation of finding sex-equal, 
egalitarian, minimum regret, or maximum cardinality stable sets for an instance of the SMP with 
unacceptability and ties. Hence, our general framework allows us to tackle a class of problems and 
requires only small adaptions to easily shift between them. 
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